Diffusive transport in spin-1 chains at high temperatures 
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We present a numerical study on the spin and thermal conductivities of the spin-1 Heisenberg 
chain in the high temperature limit, in particular of the Drude weight contribution and frequency 
dependence. We use the Exact Diagonalization and the recently developed microcanonical Lanczos 
method; it allows us a finite size scaling analysis by the study of significantly larger lattices. This 
work, pointing to a diffusive rather than ballistic behavior is discussed with respect to other recent 
theoretical and experimental studies. 
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Introduction.- Recently, numerous experiments on 
uasi-one dimensional (ID) spin-1/2 compounds 0, 0, 
"~| 0| have confirmed highly anisotropic thermal trans- 
port along the direction of the magnetic chains and a 
large contribution to the thermal conductivity due to the 
magnetic interactions. This is in agreement with early 
theoretical proposals |E °f ballistic transport in spin- 
1/2 Heisenberg antiferromagnetic chains (HAFM), that 
was recently related to the integrability of this system 
d, 0, El 113 ■ These developments promoted the the- 
oretical study of several models, as spin-1/2 frustrated 
chains, ladders and higher spin systems, using numeri- 
cal method s lip. ll3L H3 | or low energy effective theories 
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On the spin-1 compound AgVP2S6 Hcj, thermal con- 
ductivity experiments revealed anisotropic transport - 
qualitatively similar to that of spin-1/2 compounds - 
while NMR [2l| concluded to diffusive spin transport 
at high temperatures and suggested a change in behav- 
ior at low temperatures. The 5 = 1 HAFM model is 
nonintegrable and its physics characterized by a finite 
excitation gap Although there has been signifi- 

cant progress in understanding the thermodynamics of 
5 = 1 compounds, there are still open questions regard- 
ing transport. In particular, theoretical analysis based on 
a semiclassical approach of the quantum non-linear sigma 
model (NLcrM) [11 Ell - the standard low energy descrip- 
tion of 5 = 1 chains and an integrable model - concluded 
to diffusive dynamics while a Bethe ansatz method cal- 
culation [23|,|24( to ballistic transport. 

The present experimental and theoretical status opens 
two perspectives that motivate this work; first, once the 
ID magnetic transport was established as a new mode 
of thermal conduction, the ongoing synthesis and study 
of novel compounds demands the theoretical character- 
ization of conductivities - ballistic vs diffusive - in vari- 
ous spin models. Second, the conjectured connection of 



ballistic (dissipationless) transport to the integrability of 
systems requires further theoretical analysis and confir- 
mation. 

In this paper, we present a numerical analysis of 
the thermal and spin transport properties of the spin- 
1 HAFM system in an attempt to obtain a first, albeit 
for finite size lattices, exact picture of the finite temper- 
ature/frequency dynamics of this prototype model. We 
focus the analysis to high temperatures in order to min- 
imize finite size effects and draw reliable conclusions on 
the thermodynamic limit. In particular, we evaluate the 
thermal/spin Drude weights, used as the criterion of bal- 
listic or diffusive transport. Additionally, we perform 
calculations for the spin u{lj) and thermal k(u>) conduc- 
tivity spectra using the Exact Diagonalization (ED) and 
the recently developed Microcanonical Lanczos Method 
(MCLM) which allows us to obtain results for 
larger systems than hitherto accessible. The data can 
be used as a benchmark in the development of analyti- 
cal theories and in the interpretation of experiments in 
spin-1 compounds. 

Model and Method.- The Hamiltonian of the spin-1 
HAFM chain is 

L 

H=jJ2SfS l+1 , (1) 

i=i 

where Si is a spin-1 operator at site / and J the exchange 
constant. We consider periodic boundary conditions and 
set J = h = ks = 1. The spin j z and energy j E current 
operators obtained from the continuity equations for the 
local spin S z and energy H are, 

L 

f = jJ2SfSf +1 -SfSf +1 , (2) 

2 = 1 
L 



(3) 
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where P are the permutations of x, y, z. 

Within linear response theory [2(| 123, |28j the real part 
of the thermal conductivity at frequency lu and temper- 
ature T is given by 

k(lu) = 2irD th 8{u) + K reg (u>), (4) 

where the regular part of the conductivity K reg is 

K reg (u) = ^tanh(^)^ J + dte* zt ({ 3 E (t),3 E }), 

~ ' " (5) 

and the thermal Drude weight D t h is obtained from 

Dth = ^L Pn\(m\j E \n)\ 2 . (6) 

Here (3 — 1/T , z = lu + irj, p n are the Boltzmann weights 
and \n) (\e n )) the eigenstates (eigenvalues), while in Eq. 
© the symbol (} denotes a thermal average. In the (3 — » 
limit we can derive the sum-rule 

r+oo 02 
J du, K (u>) = ^-((j E ) 2 )=I, (7) 

suggesting that a measure of the ballistic contribution 
to the conductivity is given by the quantity 2irD t h/I- 
The corresponding equations for the regular part of the 
spin conductivity a reg {oj) and Drude weight D, can be 
obtained from Eqs. l">l(il) above by replacing j E by j z 
and dividing them by (3. 

Drude weight data are obtained by using ED which 
restricts us to system sizes up to L = 12 sites. We use 
the translational and spin symmetries of our system to 
perform the calculation in subspaces of momentum k and 
magnetization S% ot . We find that the results obtained in 
the k = 0, Sf 0t — subspace for L — 12 (space dimension 
~ 6500) are very close to those obtained by diagonalizing 
the entire Hilbert space. 

For the high temperature k(lu) and cr(u>) calculations 
we employ the MCLM method [Til Ejf which allows us 
to obtain results for systems up to L = 18 sites. The 
spectra calculated using this method include the Drude 
weight as a low frequency peak with width of the order of 
the frequency resolution of the method; notice however 
that this contribution is negligible for the larger systems 
we study as it follows from the finite size scaling of the 
Drude weights (see Figs. 1,4). Here we have used ~ 
1000 Lanczos steps for the first Lanczos procedure and ~ 
4000 Lanczos steps for the continuous fraction expansion 
which results in an lu resolution of ~ 0.01. 

Thermal conductivity. - In Fig. ^we show the temper- 
ature dependence of the thermal Drude weight for several 
system sizes L. D t h is vanishing at T = while at high 
temperatures it has a simple (3 2 dependence. A nonzero 
Drude weight is generally expected for systems with size 




FIG. 1: The thermal Drude weight as a function of tem- 
perature for several system sizes. In the inset, the scaling of 
Dt h /P 2 for p -» for odd and even number chains. 

less than the mean free path of the magnetic excitations. 
In the (3 — > limit, as L increases, D t h/ ' 1 decreases - 
seemingly exponentially fast - and appears to scale to 
zero in the thermodynamic limit as seen by the curves in 
the inset of Fig. (for chains with even and odd number 
of sites). Our data therefore suggest diffusive thermal 
transport for the spin-1 HAFM chain. 

We now apply the MCLM method to calculate the lu 
dependence of the thermal conductivity k(lu) in the high 
temperature limit as shown in Fig. [21 For frequencies 
lu > 0.05 the conductivity obtained for L = 18 has prac- 
tically converged to the L — > oo limit while in the low 
frequency regime there is a remaining size dependence. 
This is partly due to the variation of the Drude weight 
which contributes to the low frequency k(lu). It is worth 
noting that the statistical fluctuations in our MCLM re- 
sults are very small, even for the smallest size system 
displayed here. A comparison of ED versus MCLM re- 
sults for L = 12 (not shown for clarity) gives satisfactory 
agreement although for this size system the statistical 
fluctuations are significant. 

On the low frequency region, it is not well described by 
a Lorentzian, as predicted by the diffusion phenomenol- 
ogy j2^|, but the overall form of k(lu) is similar to that 
found in the 5=1/2 ladder model [lj] and other low 
dimensional models [30|: it suggests that this may be a 
generic behavior of conductivity spectra in such systems. 
From this curve we can also extract an estimate |l 1| of 
the high temperature k^c — k(lu —>■ 0) thermal conductiv- 
ity, Kd c ~ 16(/9J) 2 ^:, assuming typical lattice constants 
O(10A) and J - 0(1000 K). 

For comparison we note that, (i) a k^c ~ 0(1 -^e?) was 
observed at temperatures below the gap - T ~ 0.2J - in 
the compound AgVP2S6 H3] and (ii) our high temper- 
ature Kdc for the spin-1 model is an order of magnitude 
larger than that of the ladder [l4|; notice that the spin-1 
HAFM has a similar low energy excitation spectrum and 
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FIG. 2: The thermal conductivity k(lo) for (3 — ► 0. 








FIG. 3: The dependence of the Drude weight for L = 10. 

thus low temperature behavior as the spin-1/2 two-leg 
ladder for J±/J ~ 0.9. 

Spin conductivity.- We now investigate the spin trans- 
port by calculating the Drude weight D and spin conduc- 
tivity a(u>). For L ^ 12, D/(3 appears to be equal to zero 
(up to numerical precision) at all temperatures. On this 
issue it is important to point out that, for faster conver- 
gence, we consider only the S z — subsector that is the 
dominant one in the thermodynamic limit. In order to 
explore the robustness of this result we apply the canon- 
ical transformation 5 ; + — > S^e 1 ^ 1 on H and j z (periodic 
in <f) with period 2it/L); the results for D as a function 
of 4>, are shown in Fig. 03 for L = 10. 

We find that D is finite for all 4>, except for cf> = and 
(\> = ir/L where it develops a sharp minimum. Curves for 
different L's show very similar </> dependence, but with 
the local minimum at cj> = ir/L becoming sharper with 
increasing L. We therefore conclude that the vanishing 
Drude weight at <f> = 0, even for small L, is an artifact 
of the periodic boundary conditions; notice that a van- 
ishing D / f3 and a nontrivial cj> dependence is also found 
in the 5 = 1/2 isotropic model. In Fig. 0|we show that 
D/0, close to its maximum value at <f> = tt/2L, is finite 
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FIG. 4: Drude weight D as a function of temperature for 
several system sizes L. In the inset, the scaling of D//3 for 
/3-»0. 
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FIG. 5: The spin conductivity cr(ui) for ft 0. 



throughout the temperature range and as shown in the 
inset, it scales to zero exponentially fast with L in the 
(3 — > limit. In contrast to the thermal Drude weight, 
D/P goes to a finite value at very low temperatures that 
can be understood considering that a <fi implies a ground 
state carrying a nonzero spin current. 

Finally, we present in Fig. 0<t(<j) evaluated using the 
MCLM method. We see that there are some statistical 
fluctuations in the data for the smaller systems while 
those for the larger systems are very smooth. The curves 
seem converged to their L — > oo limit for u> > 0.2. The 
main characteristics of our results is the appearance of a 
local maximum at w ~ 1 and a minimum at to < 0.2. The 
later disappears with increasing system size while again 
we see no signs of a Drude peak in the ct(uj) curves. It is 
interesting to note that the local maximum feature has 
also previousl y b een observed in a study of correlations 
of the NLcrM [13. 

Discussion.- The overall picture emerging from the 
presented numerical data shows that the high temper- 
ature spin and energy transport of the spin-1 HAFM 
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chain is characterized by finite dc values, vanishing Drude 
weights, a smooth frequency dependence (though not of a 
Lorentzian form) and thus non-ballistic character. This 
behavior is compatible with the assumption of normal 
transport in nonintegrable models, it is qualitatively sim- 
ilar to that of spin- 1/2 antiferromagnetic ladder and in 
contrast to the ballistic transport of the integrable spin- 
1/2 version. On this point we should mention that in 
the isotropic spin-1/2 model D also seems to vanish ||; 
however o~(uj — > 0) might diverge |2fij and in any case, in 
the easy-plane anisotropic 5=1/2 model D is finite, (in 
contrast to preliminary results on the anisotropic 5 = 1 
model). On the other hand, for 5 — 1/2, D t h is clearly 
finite, as the energy current operator commutes with the 
Hamiltonian, again in contrast to the 5 = 1 case. 

From our data in Figs. [2Jand[S]we can also extract the 
spin T> s = adc/x ~ l-4/3/|/3 ~ 2.1 and thermal T> t h = 
Kdc/C ~ 7.5/3 2 /|/3 2 ~ 5.6 diffusion constant (in units 
of J/K), where x is the static susceptibility and C the 
specific heat. In comparison, a standard (3 — > moment 
analysis @] gives, V s = ^/2?r5(5 + l)/3 - 2.1 and V th = 
v/tt5(5+ l)/3/(l - 3/45(5 + 1)) - 2.3; the agreement 
for T> s is excellent (also probably fortuitous considering 
the quantum character of the 5 = 1 system) while an 
enhanced value is found for T>th- 

Regarding the low temperature behavior, the limita- 
tion of our calculation to small size systems (thus a sparse 
low energy spectrum) does not allow us to make any re- 
liable statements and in particular, to study an even- 
tual change of transport from diffusive to ballistic as sug- 
gested by the experimental results for AgVP2S6 Hj. Yet, 
there exist spin-1 compounds known, with weak values 
of J, for which these data are directly relevant in the 
interpretation of transport experiments. A crucial issue 
remains however for future studies, namely the disentan- 
glement of the spin-phonon from the intrinsic spin-spin 
scattering contribution to diffusion. 

Finally, on the low energy NLcrM approach [IH.l23ll24]| . 
this high temperature study cannot shed light on the 
issue of ballistic vs. diffusive behavior. If it is concluded 
that the NLcrM predicts diffusive transport then there is 
continuity with the present j3 — ► data. If, on the other 
hand, ballistic transport (perhaps due to the integrability 
of the NLcrM) is found then, the omitted "irrelevant" 
terms (for the thermodynamics) could result to a diffusive 
behavior at all temperature scales. 
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